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Neutron stars are not only described by their mass and radius, but also by how fast they spin 
(moment of inertia) and how much they can be deformed (Love number and quadrupole moment). 
These depend sensitively on the star's internal structure. We find universal relations between the 
moment of inertia, the Love number and the quadrupole moment that are independent of the equa- 
tion of state. These relations can be used to learn observationally about neutron star deformations, 
break degeneracies in gravitational wave observations, and test General Relativity independently of 
nuclear-structure. 
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Introduction. One of largest uncertainties in nuclear 
physics is the relation between energy density and pres- 
sure at incredibly high densities, the so-called equation of 
state (EoS) of nuclear matter. The interior structure of 
neutron stars (NSs) depends sensitively on its EoS, and so 
does its exterior properties, such as its mass and radius, 
as well as how fast the star can spin, characterized by 
its moment of inertia, and how much the star can be de- 
formed, characterized by its quadrupole moment and its 
tidal Love number. The astrophysical observation of NSs 
and their exterior properties may allow us to infer the NS 
EoS [lM|. For example, the observation of X-ray bursters 
and low-mass X-ray binaries has allowed for the .simulta- 
neous determination of the mass and radius to O(10)% 
accuracy Observations of millisecond pulsars, such 

as the double binary J0737-3039 @, 01: ma-Y allow for the 
measurement of the moment-of-inertia to the same accu- 
racy d, II] . Gravitational wave (GW) observations from 
binary NS inspirals with second-generation ground-based 
detectors, such as Adv. LIGO, Adv. VIRGO and KA- 
GRA, may allow for the measurement of the tidal Love 
number jl0l - ll5l |. 

None of these observations, however, are currently ac- 
curate enough to select between the many different EoSs 
that have been proposed, which then leads to degenera- 
cies in the extraction of NS information from new obser- 
vations. For example, GW observations cannot extract 
the individual NS spins, because these are degenerate 
with the NS quadrupole moment. Similarly, GWs from 
NS binary inspirals cannot be easily used to test General 
Relativity (GR), again due to degeneracies in the equa- 
tion of state |16l - [l8| . We here find a way to uniquely 
break these degeneracies through universal I- Love- Q re- 
lations, ie. analytic relations between the moment-of- 
inertia /, the tidal Love number A^*"^' and the quadrupole 
moment Q that are essentially insensitive to the NS EoS. 

The reason why these universal relations hold is two- 
fold. On the one hand, we find a strong evidence that 
these quantities depend most sensitively on the NS inte- 
rior structure close to the crust, precisely where our EoS 
ignorance is minimal and realistic EoSs agree. On the 
other hand, the efFacement principle and no hair theo- 
rem hold in GR, which state that as a body approaches 



the test particle limit, the motion and the exterior field 
of a body does not depend on its internal structure. 

These universal relations have applications in several 
different fields. On an observational astrophysics front, 
the measurement of a single member of the I-Love-Q trio 
would automatically provide information about the other 
two, even when the latter may not be accessible to ob- 
servation. Such a measurement would provide all the 
information necessary to describe the exterior proper- 
ties of tidally-deformed and slowly-rotating NSs at lin- 
ear and quadratic order in spin. On a GW front, the 
LLove-Q relations would break the degeneracy between 
the quadrupole moment and the spins in GWs emitted 
during binary NS inspirals. Given a GW detection from 
such a source with a second-generation ground based de- 
tector, such as Adv. LIGO, the GW data analysis com- 
munity may then be able to measure the individual NS 
spins to about 0.01 in dimensionless units. On a funda- 
mental physics front, the TLove-Q relations will allow, 
for the first time, for tests of GR in the NS strong-field 
that are not only theory-independent, but also EoS in- 
dependent. In particular, independent measurements of 
two members of the I-Love-Q trio would constrain vi- 
olations of the strong-equivalence principle, one of the 
pillars of GR, orders of magnitude more effectively than 
Solar System constraints. 

Universal Relations. Consider an isolated, slowly- 
rotating NS that is described by its mass M, the mag- 
nitude of its spin angular momentum J and angular ve- 
locity r2, its (spin-induced) quadrupole moment Q and 
its moment of inertia / = J/51. Let us introduce di- 
mensionless quantities / = I/IvP and Q = —Q/{M^\^), 
where x = J/M'^ is the dimensionless spin parameter. 
Physically, / determines how fast a body can spin, given 
a fixed J, while Q encodes how much a body can be 
quadrupolarly deformed. These quantities are deter- 
mined by solving the perturbed Einstein equations in the 
slow-rotation limit to first and second order in spin, re- 
spectively [H, H^, given an EoS. 

In the presence of a companion, a NS will also be 
quadrupolarly deformed. The quadrupole moment ten- 
sor Qij determines the magnitude of this deformation 



and it can be written as Qi 



-AC'^^^fy, where A^"'^) is 
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FIG. 1: (Top Left and Right) Universal I-Love and Q-Love relations for various EoSs, together with fitting curves (solid curves). 
On the top axis, we show the corresponding NS mass M with an APR EoS. (Bottom Left and Right) Fractional errors between 
the fitting curve and numerical results. 



the tidal Love number and £ij is the quadrupole (grav- 
itoelectric) tidal tensor that characterizes the source of 
the perturbation [lO, 14, SOj, [2l|. Let us introduce the di- 
mensionless tidal Love number A''"*-' = A^*'''^-^^; which 
physically characterizes the tidal deformability of a NS in 
the presence of its companion's tidal field. A''"^-' can also 
be calculated by treating the tidal effect of the compan- 
ion star as the perturbation to the isolated (non-rotating) 
NS solution 

We here derive universal relations between /, Q 
and that are essentially insensitive to the EoS. 

We consider 4 different realistic EoSs: APR j2^ . 
SLy (23t . Lattimer-Swesty with nuclear incompressibil- 
ity of 220MeV (LS220) ^ and Shen ^^i- For 
the latter two, we adopt a temperature of O.OlMeV and 
an electron fraction of 30%. We also consider a simple 
n = 1 polytropic EoS, with p = Kp^^^/^. We assume 
the NS is uniformly and slowly-rotating, with isotropic 
pressure. 

The top panels of Fig. 1 present the LLove and Q-Love 
relations for various EoSs. Observe that these relations 
hold universally, essentially independently of the EoSs. 
Such relations can be fitted with a polynomial on a log- 
log scale [l^l, shown here with solid black curves. The 
bottom panels of this figure show the fractional errors 
between the fitted curves and the numerical results. A 
similar universal relation holds between I and Q [20j . 

The LLove-Q relations hold mainly for two reasons. 
One of them is that these quantities depend the most 
on the NS internal structure only near their crust, where 
our ignorance of nuclear physics is minimal and realis- 
tic EoSs agree. We verified this explicitly by computing 
the LLove-Q relations for NSs with n = 2, 2.5 and 3 
polytropic EoSs, finding that indeed the I-Love-Q curves 
begin to deviate away from those shown in Fig. 1 as one 
increases n. The other reason is related to the effacement 



principle [H, [13, HI] , a consequence of the no- hair and 
strong equivalence principles. This principle says that 
the motion of a body and its gravitational field depend 
less and less on its internal structure as the compactness 
of the body, the ratio of its mass to its radius, increases. 
For NSs, the compactness is large enough (of order a 
tenth) that any possible EoS variability in the I-Love-Q 
relations mostly effaces away. We verified this explic- 
itly by computing the I-Love-Q relations for NSs in a 
modified gravity theory where the effacement principle 
is violated, finding that indeed the universal relations do 
not hold then. 

Application to Observational Astrophysics. Millisecond 
binary pulsars have the potential to measure / with 10% 
accuracy in the near future Q . The moment of inertia 
may be measurable because it induces additional preces- 
sion of the periastron, as well as precession of the angular 
momentum vector and the NS spin vectors. 

Given an observational measurement of the moment of 
inertia to 10%, the I-Love-Q relations automatically pro- 
vide a secondary measurement of the tidal Love number 
and of the quadrupole moment to roughly the same ac- 
curacy. Notice, however, that the tidal Love number and 
the quadrupole moment would not be easily directly ob- 
servable with millisecond binary pulsars. Although they 
do induce additional precession, their effect is suppressed 
relative to that of the moment of inertia, by various pow- 
ers of the binary's orbital velocity to the speed of light. 

Application to GW Astrophysics. Interferometric GW 
detectors are most sensitive to the GW phase of the sig- 
nal. For waves emitted during NS binary inspirals, the 
GW phase contains a term proportional to the NS's spin- 
induced quadrupole moment, Qi and Q2, and another 
term proportional to its tidally-induced quadrupolar de- 
formation A^*"^^ and Aj*"^^. The former enters with a 
factor proportional to v'^/c'^ [2^, while the latter en- 



3 



ters with a factor proportional to v^^ /c^'^ [l^ relative 
to the leading-order term. Since by Kepler's third law, 
the orbital velocity is related to the GW frequency via 
V oc /^/"^j each term has a distinct frequency dependence 
that in principle makes them non-degenerate. 

The NS quadrupole moment is degenerate with the 
NS's individual spins, because there is a spin-spin inter- 
action term in the GW phase that enters at the same 
order in v/c as the quadrupole one [1^ [s^l- Such a de- 
generacy prevents us from independently extracting the 
quadrupole moment and the individual spins from a GW 
detection with any second- or third-generation GW de- 
tector. The Q-Love relation, however, can be used to 
break this degeneracy, by rewriting Q as a function of 
^(tid)_ j£ ^YiQ Love number can be measured with a GW 
detection, then one can also separately measure the spins, 
which would be practically impossible without the Q- 
Love relation. 

Let us then consider the accuracy to which NS pa- 
rameters can be measured, given a binary NS obser- 
vation. We estimate this accuracy by carrying out a 
Fisher analysis with two different post-Newtonian NS 
waveforms [20| that differ only in the choice of system 
parameter. In one case, we choose the parameters [T5| 

{9\} = i\nMMv,/3,DL,t,,(j,„~xf''^), where M = 
(7711+7712)77^^^ is the chirp mass, rj ~ 7^170,2/ (70,1+7712)^ is 
the symmetric mass ratio, 7711 and 777,2 are the component 
masses, Z?l is the luminosity distance and and (/)c are 
the time and phase of coalescence. The quantity /? is a 
certain combination of the individual NS spins [30| , while 

Ai"^' is the averaged, dimensionless tidal Love number. 
In the other case, we choose the parameters [1^ {^'bI ~ 
(lnX,(5m,Xs,Xa,-DL,<c,'/>c,0s(Ai*"^^), aI*"^-*), where Xs 
is the averaged, dimensionless NS spin, Xa is the dif- 
ference of the dimensionless spins divided by two, Qs is 
the averaged quadrupole moment and we have used the 
Q-Love relation to express Qs as a function of aI*"*' . 

Figure 2 shows the measurement accuracy of different 
parameters, given a binary NS observation at Dj^ = 100 
Mpc with Adv. LIGO and a signal-to-noise ratio of 
SNR K, 30 for three different systems: (i) {mi,m2) = 
(1.45,1.35)Mo, XI =X2, (ii) (mi, "12) = (1.45, 1.35)Mq, 
Xi = 2x2 and (iii) {mi, 1^2) = (1.4, 1.35)M0, xi = X2- 
The solid red line corresponds to the measurement ac- 
curacy of P with waveforms parameterized by Oai while 
the other lines correspond to the measurement accuracy 
of Xs and Xa with waveforms parameterized by Ob- Ob- 
serve that A/3 and Axa are almost identical for all sys- 
tems, since they are dominated by their priors. Observe 
also that these errors are so large that the relative frac- 
tional errors can be larger than 100%. On the other hand, 
thanks to the Q-Love relation with the parameterization 
of 9b, one can determine the averaged spin Xs to ^ 0.01. 

Application to Fundamental Physics. NS observations 
allow for tests of GR [l^ when the gravitational field is 
much stronger than in the Solar System, where Einstein's 
theory has already been stringently tested [l^ [13 ■ Un- 
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FIG. 2: Measurement accuracy of spin parameters /9, Xs and 
Xa with Adv. LIGO given a detection at a luminosity distance 
of 100 Mpc with SNR ^ 30. We consider three different NS 
binaries, labeled by (i), (ii) and (iii), as described in the text. 
A/3 and Axa are almost identical among all systems, since 
they are dominated by their priors (|/3| < 0.2 and \xa\ < 0.1). 



fortunately, however, these tests arc most of the time not 
effective because of degeneracies between modified grav- 
ity effects and the EoS. The I-Love-Q relations can be 
used to break this degeneracy and thus allow for model- 
independent and EoS-independent tests of GR with NS 
observations, which could be thought of as null tests of 
the strong-equivalence principle. 

The most robust GR test would require at least two 
independent measurements of any two quantities in the 
I-Love-Q trio. Given a single measurement of any one 
of them, the I-Lovc-Q relations give us what the other 
two values must be in GR. The second measurement 
can then be used as a redundancy test. A disagree- 
ment between this measurement and the I-Love-Q pre- 
dicted value would indicate a model-independent and 
EoS-independent GR deviation. 

As an example, let us assume that a millisecond bi- 
nary pulsar observation has led to a measurement of / to 
10% accuracy [1,0, while a GW observation has led to a 
measurement of A*-'"^^ to roughly 60% for a system with 
similar masses [lol - [l5l [20j . Such a measurement would 
lead to an error box in the I-Love plane, centered about 
the measured value, as shown in the left panel of Fig. 3. 
Any modified gravity theory must then be such that its 
I-Love relation runs through this error box. This panel 
shows not only the I-Love relation in GR, but also this 
relation in a certain modified _gravity theory (dynamical 
Chern-Simons (CS) gravity [3l|) for a variety of EoSs. 
The requirement that the CS I-Love relation be inside of 
the error box automatically places a stringent constraint 
on this theory, which would be approximately 6 orders of 
magnitude stronger [20j than current Solar System con- 
straints [33 |. 

Even if one were not able to measure the moment of in- 
ertia in the near future, one can still test GR through the 
relation between the tidal Love number and the NS com- 
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FIG. 3: (Left) Possible error box in the I-Love plane, given two independent observations of the moment of inertia and the tidal 
Love number, shown with a black cross. The black solid line shows the LLove relation in GR, while all the other lines show the 
same relation in dynamical CS gravity. (Right) Possible error box in the Love-C plane, given two independent observations of 
the Love number and the NS compactness, shown with a black cross. The different curves show the Love-C relation in GR for 
several EoSs. 



pactness C. Let us then assume that a GW observation 
has measured A*^*"^^ to roughly 70% [1^ and that future 
low-mass X-ray binary observations have measured C to 
5% for a system with similar masses 0, [^. Such obser- 
vations would lead to an error box in the Love-C plane, 
centered about the observed value, as shown in the right 
panel of Fig. 3. The Love-C relation, however, is depen- 
dent on the EoS, as shown by the different EoS curves 
in this figure. However, the vertical distance between all 
curves is much smaller than the error in the measurement 
of the tidal Love number, making the Love-C relation 
effectively EoS-independent. The requirement that any 
Love-C relation run through this error box allows for a 
model-independent and an effectively EoS independent 
test of GR. 

Discussions. The LLove-Q relations open the door to 
exciting applications in astrophysics, GW theory and 
fundamental physics. We have here performed a cur- 
sory study of possible applications, but these could be 
followed up by much more detailed analysis. For exam- 
ple, the measurement accuracy of GW phase parameters 
was here estimated via a Fisher analysis, but this could 
be improved through Bayesian methods [H, [s^]- One 
could also extend these tests to systems that do not have 
the same masses. Millisecond binary pulsar observations 
may indeed not observe NS systems with the same masses 



as GW observations. We have here verified, however, 
that all the applications discussed above are robust, even 
when all the NS masses differ by about 10% [20j . 

The analysis of the LLovc-Q relations presented here 
opens up the road for multiple follow-up studies. For ex- 
ample, one could determine whether these relations hold 
for NSs with anisotropic pressur e [35 | . large internal mag- 
netic fields and rapid rotation |36l439f . One could also 
investigate whether there are universal relations between 
other quantities, such as the f- and w-modes of NS oscil- 
lations gmil. 
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